Giant magnetic anisotropy of the bulk antiferromagnets IrMn and IrMna 
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Theoretical predictions of the magnetic anisotropy of antiferromagnetic materials are demanding 
due to a lack of experimental techniques which are capable of a direct measurement of this quan- 
tity. At the same time it is highly significant due to the use of antiferromagnetic components in 
magneto-resistive sensor devices where the stability of the antiferromagnet is of upmost relevance. 
We perform an ab-initio study of the ordered phases of IrMn and IrMua, the most widely used 
industrial antiferromagnets. Calculating the form and the strength of the magnetic anisotropy al- 
lows the construction of an effective spin model, which is tested against experimental measurements 
regarding the magnetic ground state and the Neel temperature. Our most important result is the 
extremely strong second order anisotropy for IrMus appearing in its frustrated triangular magnetic 
ground state, a surprising fact since the ordered LI2 phase has a cubic symmetry. We explain 
this large anisotropy by the fact that cubic symmetry is locally broken for each of the three Mn 
sub- lattices. 

PACS numbers: 75.30.Gw 75.50.Ss 71. 15. Mb 71.15.Rf 



While the magnetic anisotropy (MA) of ferromagnets 
is a well investigated quantity, both experimentally as 
well as theoretically, it is much less understood in case 
of antiferromagnets. This lack of knowledge is on the 
one hand due to a lack of experimental techniques which 
are capable of a direct measurement of this quantity. On 
the other hand, theoretical first principles calculations of 
magnetic anisotropy effects are quite challenging as they 
require the use of fully relativistic spin density functional 
theory. 

Interest in the MA of antiferromagnets comes from the 
fact that these compounds arc important components of 
GMR sensors used, e.g., in hard disc read heads. An- 
tiferromagnetic materials are employed in these devices 
to form antiferromagnet/ferromagnet bilayers exhibiting 
exchange biaa^, a shift of the hysteresis loop of the fer- 
romagnet, providing a pinned layer which fixes the mag- 
netization of the reference layer of a GMR sensor. The 
stability of the antiferromagnet is most crucial for the 
stability of exchange bias and hence the functioning of 
the device^i^. Industrially the antiferromagnet IrMn is 
widely used because of the large exchange bias and ther- 
mal stability that can be obtained with this material. 

From experimental investigations of the exchange bias 
effect it is concluded that IrMn must have a rather large 
MA. Recent estimates of the MA of IrMn concerned 
the mean blocking temperature Tb, the temperature at 
which the exchange bias shift changes sign upon thermal 
activation. From T^ the intrinsic MA can be inferred if 
the particle size distribution is known; such a procedure 
has recently been reported and the room temperature 
MA energy of IrMn was estimated at 5.5x lO^erg/co^ and 
even 2.8x lO^erg/co^ depending on the seed layer and, 
consequently, the texture of the IrMn. 

In this letter, we address several features of the MA 



of IrMn alloys starting from first principles. In terms of 
simple symmetry considerations we predict the form of 
the MA energy that we fully confirm using ab-initio cal- 
culations providing also the strength of the MA, i.e., the 
relevant MA constants. To our best knowledge, for frus- 
trated antiferromagnets, such as IrMna, this is the first 
theoretical prediction of the MA in the literature. Our 
most remarkable observation is the surprisingly strong, 
second order MA of IrMna resulting from the fact that 
the cubic symmetry is locally broken for each of the three 
Mn sub-lattices. We are also able to attribute contribu- 
tions of the MAE related to on-site and two-site exchange 
anisotropy terms, a very important issue for finite tem- 
perature magnetism^ii. Such a separation is inevitably 
important for the purpose of subsequent simulations to 
study exchange-bias systems based on these compounds, 
for example in determining the scaling behavior of the 
MA energy. 

Self-consistent calculations are performed in terms of 
the fully relativistic screened Korringa-Kohn-Rostoker 
(SKKR) method^. Within this method, spin-polarization 
and relativistic effects, in particular, spin-orbit coupling 
are treated on equal theoretical footing by solving the 
Kohn-Sham-Dirac equation. The local spin-density ap- 
proximation as parametrized by Vosko et al^ was ap- 
plied; the effective potentials and fields were treated 
within the atomic sphere approximation with an angu- 
lar momentum cut-off of Imax = 2. For the LIq IrMn 
alloy we used the lattice constants a — 3.855A and c = 
3.644jy^, while for the LI2 IrMna alloy an fee lattice with 
a = 3.785A was considere d^^'^^ . For the self-consistent 
calculations we fixed the orientations of the magnetic 
moments on the Mn atoms according to the magnetic 
ground-state configurations reported previously in the 
literature, namely, a checkerboard coUinear AF structure 
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for LIq IrMniSii^ and a triangular (Tl) state within the 
fcc(lll) planes for LI2 IrMnqfii— . We obtained vanish- 
ing spin-polarization at the Ir sites, whereas spin mag- 
netic moments of 2.63 /is and 2.66 /.ts at the Mn sites 
for IrMn and IrMna, respectively. These values are in 
satisfactory agreement with earlier first principles calcu- 
lationsiSii^. 

We start our study of the MA by symmetry consider- 
ations based on the following effective spin-Hamiltonian 
(energy per unit cell), 
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where Sa is the spin-vector of the Mn sub-lattice la- 
beled by a; n = 2 for Llg IrMn and n = 3 for LI2 
IrMna. Dab are (traceless) symmetric matrices repre- 
senting anisotropic two-site (exchange) coupling and Ka 
are on-site anisotropy matricesJ^ Note that all the pa- 
rameters in Eq. ^ are defined as sums over sites in the 
sub-lattices, e.g.. Jab = J^jeb'^v for i G a (j ~ i ex- 
cluded), Jij being the isotropic intersite interactions. In 
case of LIq IrMn, tetragonal symmetry implies. 
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with Dii — D22 — D and D12 = D' . Rotating an 
antiferromagnetic configuration around the (100) axis. 
Si = (0, sin (ys, cos (ys) and S2 ~ ~Si, a simple orien- 
tation (if-) dependence of the energy can be derived, 
E{ip) = E{0) + Kcs sin^ ip, introducing an effective uni- 
axial MA constant per unit ceU, Kcs = 2K + |(£)' — D). 

In order to calculate E(ip) from first principles wc 
adopted the so-called magnetic force theorem^ in which 
the previously determined self-consistent effective po- 
tentials and fields are kept fixed and the change of 
total energy of the system with respect to tp is ap- 
proached by that of the single-particle (band-) energy. 
The values for E{ip) from these calculations could be very 
well fitted with K^ft = — 6.81meV, in very good agree- 
ment with the theoretical value reported by Umctsu et 
alJ^ and also with the easy-plane anisotropy observed 
experimentalljii^. Furthermore, the on-site MA constant 
K in Eq. ([2]) can be expressed as. 
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where i labels any Mn site (see Ref. [ij for details). 
Importantly, an overall coUinear magnetic arrangement 
along the (100) axis has to be considered in these cal- 
culations. By using the above formula we obtained 
K = —2.94 meV implying that in this system the MA 



energy is dominated by the on-site anisotropy, i.e. the 
third term in Eq. p]). 

In the case of LI2 IrMna each of the three Mn atoms 
in a unit cell exhibits a second order MA due to local 
tetragonal symmetry. However, as indicated in Fig. [T] 
with different symmetry axes that have to be accounted 
for in Eq. ([T|) by suitable transformations of the matrices 
in Eq. C3 rotational symmetry around the (111) axes 
implies Dii = £'22 = -D33 = D and D12 = D23 = i^si = 
D' . Clearly, for a ferromagnetic state of the system such 
a Hamiltonian would yield a vanishing MA energy. 




FIG. 1: (color online) Sketch of the IrMns unit cell. Dark 
spheres represent three Mn atoms corresponding to the anti- 
ferromagnetic sub-lattices. The solid arrows indicate the local 
easy axes and the dotted arrows indicate the spin direction in 
the Tl ground state. 

This second order MA becomes, however, evident if all 
the spins forming the Tl ground-state are rotated around 
the (111) axis. Straightforward calculations show that 
E{ip) follows again a sin^ (/3-dependence with an effective 
MA constant, KcG = 2K +^{D + D'). Our first princi- 
ples calculations reproduced well the proposed functional 
form oiE{tp) with a value of i^eff = 10.42 meV, see Fig. [51 
Thus we conclude that the MA constant for LI2 is almost 
twice as large in magnitude than for LIq IrMn. 

We confirm the validity of the spin-Hamiltonian Eq. ([1]) 
for L\2 IrMna by applying two additional rotations of the 
spin-system. First, we repeat the rotation around the 
(111) axis by simultaneously interchanging the orienta- 
tions of the spins at the Mn sites 2 and 3. It should be 
mentioned that this triangular spin-structure (say, T2) 
corresponds to a chirality vector. 
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that is just the opposite of the chirality vector related 
to state Tl. Note also that k is normal to the plane 
of the moments and the normal component of this vec- 
tor (chirality index) k for state Tl is k = 1, while for 
state T2 k = — 1. Whilst by considering only the first 
(isotropic) term in Eq. ([T]) the energy of the these two 
states is identical, the anisotropy terms lift this degener- 
acy. Interestingly, rotating the spins in state T2 around 
the (111) axis does not induce changes in the energy of 
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the system. This is confirmed by our calculations up to 
an absolute error of 2 ^eV. Furthermore, the energy of 
state T2 should be higher by —Kes/2 than the energy 
minimum of state Tl (ip = 0). From our calculations 
we found this difference to be 5.22 meV, fitting nearly 
perfectly to the previously determined MA constant. 

Our last test to Eq. ([T|) referred to rotating the spins in 
state Tl around the (110) axis. As compared to all the 
previous cases, this rotation implies a quite complicated 
form of E{(p), 

E{ip) = E{Q) + ^(2 + sm^ip-2cosip 
8 V 

— 2V2sin(y9(l — cos<^) ) . (5) 



In Fig. [2] we also plotted the results of this calculation 
together with the fit function as above. Reassuringly, this 
fimction describes E{ip) well for the whole range of ip with 
the MA constant as obtained before (ii'cff=10.42 meV). 
Note that for the rotation around the (110) axis a,t (p ~ 
109.47° the energy of the ground-state is regained. This, 
however, is not surprising since by this rotation we obtain 
a Tl state lying in a plane normal to the (111) direction. 
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FIG. 2: (color online) Calculated change of energy of the LI2 
IrMns system when rotating the triangular Tl spin structure 
around the (111) axis (circles) and the (110) axis (squares). 
The solid lines display appropriate fits to K^s sin^((/?) and the 
function in Eq. ([5]), respectively. 

In order to calculate the on-site anisotropy parameter, 
K, we again applied Eq. ([3]) by using a ferromagnetic 
reference state oriented along the (100) direction. Al- 
though we used the effective potentials calculated from 
the Tl ground-state, because of the large difference in 
the spin-configuration of the reference state and the true 
ground-state, we expect just a rough estimate on K . The 
obtained value, K ~ 1.06 meV, indicates at best that, 
unlike the LIq IrMn alloy, in this system the MAE is 
mainly governed by two-site anisotropy, i.e., the second 
term in Eq. P^. 

In the second part of this Letter we present results of 
finite temperature simulations on the magnetism of IrMn 



compounds. With this purpose we construct a simplified 
effective spin model based on our first principles calcula- 
tions, 



(6) 



where arc isotropic Hciscnberg exchange parame- 
ters and the effective on-site anisotropy parameters k = 
Kcff/2, merging thus the effect of two-site anisotropy 
terms. Here, fii are unit vectors along the local uniaxial 
symmetry axes: for LIq IrMn is perpendicular to the 
Ir (or Mn) planes, for LI2 liMn^ fii is different for each 
of the three Mn sub-lattices, see Fig. [T] We calculated 
the parameters, Jy, by using the widely adopted torque 
methodic as extended to rclativistic calculations^^ 
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FIG. 3: (color online) Isotropic exchange interactions, Jij, 
between the Mn atoms in IrMn alloys calculated from the 
corresponding ground state magnetic configurations by using 
the torque method.— 

For both alloys, the calculated exchange interactions 
are shown in Fig. [3] as a function of the distance be- 
tween the Mn atoms. The two sets of interactions show 
obvious similarities: large antiferromagnetic (negative) 
nearest neighbor interactions, sizable oscillating interac- 
tions up to about Rij = 6 A, while a strong damping for 
larger distances. Note that double (multiple) values for 
some Rij 's appear due to the different symmetry (neigh- 
borhood) of pairs with the given separation. In case of 
LIq IrMn these 'degeneracies' are mostly resolved via 
tetragonal distortion of the lattice. In good compari- 
son with other theoretical work s ^°d^ from the calculated 
Jij's the mean-field estimates for the Neel temperatures, 
T/v = 1398K and 1222K, can be obtained, respectively. 

The model Eq. ^ is simulated by solving the Landau- 
Lifshitz- Gilbert (LLG) equation with Langevin dynam- 
ics, calculating thermal equilibrium properties in the long 
time (and high damping) limit. The methods we use 
are described in detail in Refii^. The main quantity of 
interest is the sub-lattice staggered magnetization, M^, 
defined as 



Ms 



1 " 

-E 



ax 



ay ' az 



(7) 



4 



where Aia = X^iGa proportional to the magnetiza- 
tion of sub-lattice a, n is the number of antiferromagnetic 
sub-lattices and ( ) denotes a thermal average. 
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FIG. 4: (color online) Staggered magnetizations, Ms, as a 
function of temperature obtained using Langevin dynamics 
over 20 ps with system sizes of 24000 sites {LI2) and 70000 
sites (ilo) and using periodic boundary conditions. 

Fig. [4] shows the order parameter, Mg, versus temper- 
ature T. Despite finite size effects, T/v can be estimated 
as 1360K for LIq IrMn and 1005K for LI2 IrMng. Note 
that though the exchange parameters in both cases have 
similar values, the critical temperature in the LI2 phase 
is significantly lower. Obviously, the frustration of the 
spin-ordering in the LI2 phase leads to a reduced T^r as 
compared to the LIq phase. The simulated critical tem- 
peratures clearly improve upon the mean field estimates 
as compared with experimentally observed Neel temper- 
atures, 1145K and 960K— , respectively. 

A further analysis of the sub-lattice magnetization vec- 
tors reveals the magnetic ground state configurations. In 
the case of ilo IrMn the Mn spins align along the (IIO) 
direction appropriate with the easy plane anisotropy for 
this material. For the LI2 system, magnetic anisotropy 



included according Eq. ^ reveals that the Tl ground 
state structure is fixed to lie in one of the (111) planes, 
with each of the Mn spins directed along the correspond- 
ing (211) directions. These spin-orientations have previ- 
ously been established by neutron scattering^i; our re- 
sults for the Neel temperature and the magnetic ground- 
state structures arc in excellent agreement with experi- 
ments, underpinning the validity of our spin model de- 
rived from first principles. 

In summary, we performed an ab-initio study for the 
ordered phases of IrMn and IrMna, the most important 
industrial antiferromagnets. The calculated Heisenberg 
exchange integrals and magnetic anisotropy constants are 
used to construct an effective spin model which is simu- 
lated using the stochastic Landau-Lifshitz-Gilbert equa- 
tion. A good agreement of the calculated Neel tem- 
peratures and magnetic ground-states with experimen- 
tal results confirmed the validity of our approach. Our 
most spectacular finding is a giant second order magnetic 
anisotropy for IrMns, leading to energy barriers of the or- 
der of Kcff ~ 3 X 10* erg/cc for rotation of the Tl groimd 
state spin-structure around the (111) axis. This uniaxial 
magnetic anisotropy is understood due to the fact that 
the cubic symmetry is locally broken for each of the three 
sub-lattices of the antiferromagnetic Tl ground-state. 

The extremely high anisotropy for the LI2 phase has 
perhaps not been measured experimentally because of the 
disordered nature of this material in thin film devices, 
where deposition by sputtering causes significant loss of 
long range crystallographic order. Our results, however, 
suggest that finer control of the crystallography will allow 
the extremely large anisotropy of these materials to be 
fully exploited, allowing, for example, antiferromagnet 
film thicknesses to be reduced without loss of exchange 
bias stabilitjsi^. 
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